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ABSTRACT 

In  the  past  few  years,  we  have  demonstrated  how  the  surface  return  measured  by  the  active 
instruments  onboard  CloudSat  and  CALIPSO  could  be  used  to  retrieve  the  optical  depth  and  backscatter 
phase  function  (lidar  ratio)  of  aerosols  and  ice  clouds.  This  methodology  lead  to  the  development  of  a 
data  fusion  product  publicly  available  at  the  ICARE  archive  center  using  the  Synergized  Optical  Depth  of 
Aerosols  and  Ice  Clouds  (SODA  &  ICE)  algorithm1.  This  algorithm,  also  allowing  to  derive  ocean  surface 
wind  speed,  has  been  extended  to  include  dense  cloud  surface  return  to  analyze  aerosol  and  cloud 
properties  above  such  clouds. 

This  low  level  data  fusion  of  CALIPSO  and  CloudSat  ocean  surface  echoes  has  been  used  by  several 
researchers  to  explore  different  research  paths.  Among  them,  we  can  cite: 

•  A  new  characterization  of  the  lidar  ratio  of  cirrus  clouds2 

•  The  analysis  of  the  precipitable  water  and  development  of  a  new  Millimeter-Wave  Propagation 
Model  for  the  W-Band  observations  (EMPIRIMA3) 

•  The  analysis  of  the  lidar  ratio  of  sea-spray  aerosols4,  and  of  Aerosol  multilayer  lidar  ratio  and 
extinction5 

•  A  contribution  to  the  retrieval  of  the  subsurface  particulate  backscatter  coefficients  of 
phytoplankton  particles6 

In  this  paper,  we  present  the  main  features  of  SODA  &  ICE,  summarizing  some  of  the  results  obtained. 
Keywords:  Lidar,  Remote  Sensing,  Ocean,  Aerosols,  Clouds 


1.  INTRODUCTION 

Accurate  measurements  of  aerosol  and  clouds  optical  properties  are  required  to  precisely  quantify  their  radiative 
effect.  Aerosols  scatter  and  absorb  solar  radiation  (direct  effect,  stronger  above  highly  reflecting  surfaces  including 
clouds)  and  alter  the  lifetime  and  development  of  clouds,  which  in  turn  also  impact  scattering  and  absorption  of 
radiation  (indirect  effect). 

Among  these  optical  properties,  the  aerosol  optical  depth  (AOD)  is  a  basic  quantity  representing  the  aerosol  load  in 
the  atmosphere.  The  main  error  sources  when  determining  the  aerosol  optical  depth  from  space  are  the  accuracy  of 
the  instrument  calibration,  the  cloud  screening  algorithm  used,  the  assumption  on  aerosol  microphysical  properties 
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and  the  different  parameterization  of  the  surface7.  Using  the  ocean  surface  echo  of  CALIPSO  and  CLOUDSAT 
altogether,  we  have  developed  a  new  methodology8, 9’ 10  that  allows  to  overcome  the  last  two  issues:  CLOUDSAT 
radar  provides  a  direct  measurement  of  the  ocean  surface  reflectance  in  the  microwave  which  combined  with 
CALIPSO  lidar  attenuated  surface  reflectance  in  the  visible  provides  a  direct  retrieval  of  the  total  column  aerosol 
attenuation,  without  any  assumption  on  their  type  or  microphysical  properties.  It  is  by  definition  collocated  with 
CALIPSO/CLOUDSAT  vertical  information  which  can  be  an  advantage  with  respect  to  passive  measurements  alone 
to  improve  cloud  screening.  Finally,  the  measurements  are  available  day  and  night  and  above  sun  glint  dominated 
region.  All  those  characteristics  make  our  methodology  an  interesting  approach  to  better  understand  the  aerosol 
radiative  effect  and  its  trends.  Therefore,  aiming  at  a  diffusion  to  the  scientific  community,  a  first  version  of  the 
product  (Synergized  Optical  Depth  of  Aerosols  and  ICE  clouds  -  SODA  &  ICE)  has  been  created  at  the  French 
Thematic  Center  ICARE. 

On  the  basis  of  recent  advancements  linked  to  the  use  of  SODA  in  ice  clouds2,  W-band  radar  water  vapor 
absorption3  and  feedback  from  the  scientific  community,  it  appears  that  the  first  version  did  not  express  the  whole 
potential  of  what  could  be  done  with  this  methodology.  This  leads  to  the  development  of  an  upgraded  version  which 
is  now  available  on  the  ICARE  website.  The  main  objective  of  this  paper  is  to  describe  the  product  and  its 
applications. 


2.  ALGORITHM  DESCRIPTION  SODA  VERSION  1 
GENERAL  OVERVIEW 

The  main  goal  of  the  SODA  algorithm  is  to  retrieve  the  optical  properties  of  optically  thin  atmospheric  features  over 
optically  dense  features.  SODA  version  1  data  product  contained  the  optical  depth  of  semi-transparent  aerosols  and 
clouds  over  the  ocean.  SODA  version  2  now  provides  the  optical  depth,  the  color  ratio,  and  lidar  ratio  of  aerosols 
and  clouds  over  the  ocean  and  over  liquid  water  clouds. 

The  concept  can  be  further  extended  over  land11  or  to  retrieve  underwater  biological  information6.  These 
applications  are  out  of  the  scope  of  this  paper  and  oceanic  applications  will  only  be  briefly  discussed  in  the 
perspectives  section  of  this  manuscript. 

SODA  version  1  is  based  on  a  previous  publication10  and  the  details  of  the  algorithm  are  presented  on  the  associated 
website1.  To  present  the  reader  with  a  better  understanding  of  SODA  version  2,  we  will  here  summarize  and 
illustrate  the  most  relevant  elements. 

The  semi-transparent  aerosol  and  cloud  optical  depth  is  retrieved  at  lidar  wavelength  through  three  steps.  Only 
profiles  over  the  ocean  are  kept  in  the  analysis.  In  order  to  emphasize  cloud  clearing,  data  profiles  are  analyzed  only 
if  the  maximum  of  the  lidar  signal  is  situated  at  ocean  level.  This  removes  profiles  with  low  level  dense  clouds  as 
well  as  profiles  totally  attenuated.  As  we  can  see  in  Fig.  1,  most  occurrences  of  the  maximum  of  the  surface  echo  for 
each  orbit  are  statistically  at  the  ocean  surface  level.  This  allows  to  find  the  peak  of  the  surface  echo. 


Altitude  (km)  Surface  signal  (km'Vsr"1) 

Fig.  1.  Left:  Distribution  of  the  altitude  of  the  maximum  of  the  surface  signal  for  one  orbit  of  nighttime  data 
(02/09/2009  around  06:53Z).  Right:  Surface  echo  used  in  the  integration  of  the  surface  signal  for  one  shot  around  the 
position  20.22N,  -78.30W. 
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After  the  position  of  the  maximum  has  been  determined,  the  integration  of  the  surface  signal  is  performed  around 
this  position  (Fig.  1).  The  upper  limit  includes  2  bins  above  the  peak  and  the  lower  limit  is  determined  by  the 
position  of  the  bin  number  572  (-0.28  km).  This  corresponds  to  the  whole  energy  contained  in  the  surface  signal, 
which  is  physically  consistent  with  the  Bidirectional  Reflectance  Distribution  Function  (BRDF)12. 

As  we  are  studying  dense  surface  signal,  saturation  may  happen  for  the  strongest  values  of  the  ocean  surface  signal8, 
10  and  has  to  be  taken  into  account.  The  way  that  has  been  chosen  to  filter  out  saturated  signals  in  SODA  VI  is  to 
exclude  surface  signal  which  maximum  exceeds  2  km_1.sr_1.  The  reason  behind  this  choice  is  the  observation  of  the 
maximum  signal  surface  distribution  (Fig.  2)  show  a  well-defined  distribution  of  values  over  this  threshold.  This 
threshold  is  higher  than  the  one  we  discussed  in  previous  studies8, 10  and  is  expected  to  correspond  to  the  presence  of 
ice  (or  other  solid  land  surface).  When  ice  is  excluded,  a  partial  saturation  is  still  possible  and  will  manifest  itself  as 
a  loss  of  linearity  of  the  surface  signal  at  low  optical  depth.  This  error  source  is  compensated  for  in  the  different 
processes  that  we  describe  below  and  included  in  the  error  budget  we  previously  determined2 


Surface  signal  maximum  (km^.sr-1) 


Fig.  2.  Distribution  of  the  maximum  of  the  surface  signal  for  one  orbit  (02/09/2009  around  06:53Z)  showing  the  main 
surface  reflectance  distribution  over  open  water  and  solid  surface. 


Monthly  global  calibration  (nighttime) 

This  procedure  allows  to  calibrate  the  method  and  through  a  lidar-radar  cross  calibration,  to  correct  any  monthly 
calibration  bias  coming  from  the  lidar  or  radar.  It  is  based  only  on  the  clear-sky  observations  and  does  not  rely  on 
any  underlying  assumption  of  surface  roughness  behavior. 

The  first  module  calculates  the  4th  order  polynomial  coefficients  linking  radar  normalized  surface  scattering  cross 
section  g’Sr  and  the  lidar  ocean  integrated  attenuated  surface  backscatter  coefficient  ysL,att  (see  reference  [10],  Fig. 
1). 

This  relationship  is  determined  on  a  monthly  basis  from  collocated  nighttime  data  and  only  for  the  clearest  lidar 
profiles  (identified  as  those  for  which  the  atmospheric  integrated  attenuated  backscatter  coefficient  -  IAB  <  0.01  sr" 
').  The  radar  data  are  corrected  from  gaseous  water  vapor  attenuation  and  from  ocean  refractive  index  p0R  variations 
with  temperature13,  10.  Water  vapor  gaseous  attenuation  is  retrieved  from  AMSR-E  integrated  water  vapor  content 
using  the  expression  derived  by  previous  research14.  The  lidar  data  are  corrected  from  gaseous  attenuation,  namely 
dry  air  extinction  and  ozone  absorption. 


Latitudinal  day  and  night  calibration 

Potential  local  calibration  biases  are  then  determined  by  calculating  the  statistical  deviations  from  the  polynomial 
relationship  previously  found  as  a  function  of  latitude.  Clear  air  profiles  (IAB  <  0.01  sr"1)  are  used  for  both  day  and 
night  observations  (see  reference  [10],  Fig.  4). 
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Fig.  3.  Atmospheric  IAB  of  all  atmospheric  features  for  one  nighttime  orbit  (02/09/2009  around  06:53Z).  The  black 
solid  line  shows  the  clear  air  threshold  (0.01  sr-1)  used  in  SODA  as  a  limit  for  clear-aircalibration. 


Optical  Depth  retrieval 

The  remaining  deviations  from  the  relationship  determined  in  the  previous  two  steps  are  then  calculated  for  each 
shot.  The  value  of  the  background  aerosol  optical  depth  in  clearest  conditions  based  on  reference  [15]  gives  a 
baseline  value  of  0.06  at  532  nm,  and  0.04  at  1064  nm  to  correct  the  residual  AOD.  We  will  discuss  this  further  with 
the  1064  nm  data.  The  corresponding  optical  depth  of  the  atmospheric  feature  is  stored  into  a  Hierarchical  Data 
Format  (hdf)  files  along  with  several  control  parameters.  The  same  naming  conventions  than  CALIPSO  files  have 
been  used  (an  example  of  file  name  is  SODA_AOD_20 10-08-01  TO  l-01-00ZN_V2-02.hdf) 

A  distinction  has  to  be  made  at  this  level  between  SODA  version  1  and  SODA  version  2.  In  version  1,  only  the 
condition  IAB  <  0.01  sr"1  was  removed  during  the  optical  depth  retrieval  process.  However,  the  other  filters 
(saturation  and  level  of  the  maximum  signal)  were  still  active  and  the  corresponding  profiles  were  excluded  from  the 
data  product.  A  different  approach  has  been  used  in  SODA  version  2,  as  we  will  now  discuss. 


3.  FEATURES  OF  SODA  VERSION  2 


OVERVIEW 

The  choice  of  removing  dense  features  in  SODA  version  1  was  aimed  at  improving  cloud  clearing.  However,  as  the 
underlying  principle  of  the  methodology  is  valid  for  ice  clouds2  in  SODA  version  2,  a  different  approach  has  been 
used. 

The  stages  of  monthly  global  calibration  and  Latitudinal  day/night  calibration  are  largely  unchanged  except  for  the 
water  vapor  correction  and  for  the  addition  of  the  calibration  of  methodology  over  liquid  water  clouds  (see  further). 
The  new  water  vapor  correction  has  been  applied  to  all  steps  of  the  algorithm  (see  reference  [3]  for  the  details  of  the 
new  parameterization). 

For  the  optical  depth  retrievals,  the  filters  have  been  removed  a  replaced  with  improved  quality  flags  and  scene 
flags.  In  version  1,  the  scene  flags  were  quite  basic  and  allowed  some  distinction  between  the  aerosols  or  clouds.  In 
version  2,  the  scene  flag  is  based  on  CALIOP  vertical  feature  mask  (VFM)  and  allows  the  user  to  select  the  scene  he 
is  interested  in  (Table  1). 

The  list  of  SODA  v2  main  changes  can  be  summarized  as 

•  New  water  vapor  correction 

•  more  evolved  scene  flag  based  on  VFM  which  provides  a  guidance  to  discriminate  aerosols  and  their  types, 
ice  clouds,  water  clouds,  etc 

•  More  quality  flags 

•  includes  dense  clouds  and  aerosols  features 

•  includes  lidar  ratio  retrieval 
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•  Includes  optical  depth  and  lidar  ratio  retrieval  over  liquid  water  clouds 
Parameters  are  determined  as  in  VI  at  532  and  1064  nm,  so  that  color  ratios  can  be  further  derived. 


Table.  1  Scene  flag  in  SODA  V 2.  Bits  are  numbered  from  right  to  left  using  a  little  endian  coding. 


Bit# 

6 

5 

4 

3 

2 

1 

0 

value 

0/1 

0/1 

0/1 

0/1 

0/1 

0/1 

0/1 

Criteria 

If  ice  detected  in 
profile 

Signification 

type  of  the  highest 
layer  :  0  = 
Randomly 
Oriented  Ice  ;  1  = 
Horizontally 
Oriented  Ice 

Stratospheric 
feature 
detected  in 
profile 

Liquid 

water 

detected  in 
profile 

Ice  detected 
in  profile 

Aerosols 
detected  in 
profile 

Clear 

profile 

Scene 

Validity 

Bit# 

14 

13 

12 

9-11 

7-8 

value 

0/1 

0/1 

0/1 

0/7 

0/3 

Criteria 

If  aerosol 
detected 

Signification 

AMSR  liquid 
water  path  >  0 
and  and  <200 

Retrievals  over 
liquid  cloud 
available 

Retrievals 

over  ocean 

available 

Subtype  of  the 
highest  layer 

Ocean  surface  type  using  the 
Lidar  surface  depolarization31: 
0=Not  determined;  1=  liquid 
water  ;  2  =  mixed  ;  3  =  ice 

Bit# 

19-21 

18 

17 

16 

15 

value 

0/7 

0/1 

0/1 

0/1 

0/1 

Criteria 

If  a  lower 
aerosol  layer 
of  a  different 
type  has 
been 
detected 

If  ice  water 
and  aerosol 
detected  in 
profile 

If  liquid  water 
detected  in  profile 

If  liquid  water 
detected  in  profile 

If  liquid  water  detected 
in  profile 

Signification 

Subtype  of 
the  lower 
layer 

All  ice  is 
above  the 
aerosols 

Aerosols  are  present 
under  the  lowest 
liquid  bin  and  no  ice 
clouds  found  under 

All  aerosols  or 
cirrus  are  above  the 
highest  liquid  bin 

All  liquid  water  bins  are 
under  the  altitude  of  the 

0  Celsius  degree 
isotherm 

The  approach  used  in  SODA  version  2  provides  the  most  flexibility  to  the  user.  It  is  however  at  the  cost  of  a  higher 
complexity  as  we  can  see  from  the  numerous  options  offered  by  Table  1.  As  not  all  users  want  to  deal  with  a 
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complex  dataset,  a  user  friendly  pre-filtered  data  product  is  planned  to  be  made  available  to  the  scientific 
community. 


OPTICAL  DEPTH  RETRIEVAL  OVER  LIQUID  WATER  CLOUDS 

Over  dense  liquid  water  clouds,  the  cross  and  co-polarized  lidar  channels  are  used  as  a  reference  to  determine  the 
atmospheric  transmission.  The  algorithm  is  here  also  based  on  the  determination  of  the  integrated  backscatter  but 
this  time  dense  clouds  are  used  as  a  target  and  lidar-only  signals  are  analysed.  As  we  focus  on  dense  clouds  as 
viewed  by  lidar  (typically  with  OD  >  3),  transmission  squared  is  smaller  than  0.005,  and  the  integrated  attenuated 

backscatter  coefficient  is  equal  to  y  =  l/2Sr|,  where  S  is  the  lidar  ratio  for  waterclouds  and  r|  the  multiple  scattering 
coefficient.  This  quantity,  determined  by  the  size  and  index  of  water  droplets,  has  been  found  to  be  constant  to  better 
than  +/-  5%  in  most  of  the  clouds. 


Depolarization  of  the  lidar  signal  into  liquid  water  clouds  to  determine  the  multiple  scattering  factor  r|  on  which  the 
integrated  cloud  signal  y  depends16.  It  had  been  adapted  in  SODA  algorithm  in  order  to  account  for  the  different 
issues  that  affect  the  signal.  More  specifically,  the  reason  behind  this  added  complexity  (described  below)  is  that  the 
data  show  a  deviation  from  theory  (1  to  1  line  on  Fig.  4  left).  Possible  reasons  are 

•  Detector  transient  response15 

•  Polarization  calibration  issues  which  could  come  from  either  the  day/night  calibration  transfer  or  the  low 
gain/high  gain  merging  process. 

The  algorithm  adjusts  the  theoretical  relationship  obtained  into  liquid  water  clouds  to  account  for  sub-micrometer  or 
non-spherical  particles16’17,18,  as  well  as  measurement  noise,  calibration  issues  and  a  non-ideal  detector  transient  (Hu 
et  al.  cloud  phase  GRL  2007),  using  a  second  order  polynomial  adjustment  as 

V  —  V  Hu(a  +  P7!  Hu)  (1) 

Where  rjHu  is  the  theoretical  relationship17,18.  To  provide  an  order  of  magnitude,  in  August  2010  the  algorithm 
determined  a  =  1.37  and  P  =  -0.69  at  532  nm  ;  a  =  1.22  and  P  =  -0.61  at  1064  nm.  This  modification  is  further 
illustrated  in  Fig.  4  (left  figure,  curved  solid  line). 


(1-5)2/(1+5)2  (Theory) 


0  50 

Latitude 

Fig.  4.  For  August  2010  nighttime  data.  Left:  observed  multiple  scattering  factor  for  dense  liquid 
water  clouds  (clear  air  above)  assuming  a  cloud  lidar  ratio  of  19  sr  as  a  function  of  the  theoretical  multiple  scattering 
coefficient  (derived  from  polarized  observations).  The  dotted  line  is  the  curve  predicted  by  theory  and  the  solid  line  is 
the  adjustment  performed  by  the  algorithm.  Right:  derivation  of  the  cloud  lidar  ratio  as  a  function  of  latitude,  this 
time  using  the  multiple  scattering  coefficient  from  Eq.  1.  The  black  dots  show  what  is  used  by  the  algorithm.  The  red 
line  is  at  19  sr. 


We  determined  the  depolarization-multiple  scattering  relationship  for  all  clouds  dense  enough  to  be  detected  on  the 
333  m  cloud  level  2  product  (which  provides  top  and  bottom  boundaries).  Clouds  where  the  maximum  of  the 
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backscatter  coefficient  in  the  profile  is  lower  than  the  cloud  base  and  the  cloud  liquid  water  path  of  AMSR-E  is  not 
higher  than  0.1  mm  were  excluded.  To  remove  the  uncertainty  on  the  cloud  water  phase,  the  cloud  top  of  the  cloud 
layer  has  to  be  situated  below  the  isotherm  0°C,  as  given  by  the  GMAO  meteorological  data  used  in  the  CALIPSO 
analysis.  A  lower  limit  in  the  integrated  backscatter  coefficient  has  been  applied  so  to  exclude  cloud  layers  which 
IAB  from  cloud  top  to  20  km  was  higher  than  the  molecular  IAB. 

For  this  case,  the  average  lidar  ratio  retrieved  is  19.4  sr  which  is  very  close  from  the  theoretical  value  (around  19.2 
sr  over  the  ocean17). 


Taking  into  account  these  constraints  in  the  analysis  of  the  multiple  scattering  factor  and  including  the  modifications 
with  respect  of  the  theoretical  relationship  for  both  lidar  wavelength,  the  lidar  equation  can  then  be  solved 
unambiguously  and  optical  depth  can  be  retrieved  with  no  assumptions  on  microphysical  properties 


LIDAR  RATIO  ALGORITHM 

The  lidar  ratio  used  in  SODA  V2  is  the  one  described  previously19.  It  has  been  chosen  because  of  its  fast 
computation  speed.  Note  that  there  are  other  ways  to  solve  the  Fernald  equation20,  2’  4.  SODA  version  1  opened  the 
door  to  detailed  research  on  the  lidar  ratio  of  specific  atmospheric  features  like  cirrus  clouds2  or  marine  aerosols4  as 
a  function  of  different  geophysical  parameters  (temperature,  wind  speed).  We  hope  that  the  supplemental  features  of 
version  2  will  allow  the  scientific  community  to  further  explore  these  topics. 


DISCUSSION  OF  RETRIEVAL  UNCERTAINTY 

1.  The  uncertainty  of  the  SODA  product  over  the  ocean  has  been  discussed  and  evaluated  in  several 
publications8,  10,  2’  3.  The  analysis  of  the  CALIPSO/CloudSat  ocean  surface  data  points  towards  an 
uncertainty  in  the  optical  depth  that  can  go  as  low  as  ±0.01 53. 

It  has  to  be  confirmed  observationally  that  the  algorithm  retrieval  is  as  accurate  as  this  theoretical  limit.  The 
observed  dispersion  with  respect  to  NASA  Langley  high  spectral  resolution  lidar  (HSRL)  flights  is  equal  to 
around  0.04.  This  is  consistent  with  our  experimental  study  based  on  the  comparison  of  SODA  with 
Infrared  observations2  where  we  have  determined  that  this  error  is  lower  than  0.05. 

MODIS  is  probably  the  best  reference  to  put  these  numbers  in  perspective  because  of  the  tremendous 
amount  of  work  that  has  been  performed  for  the  product  validation.  For  aerosol  optical  depth  iaerL ,  MODIS 
uncertainty  over  the  ocean  is  At aerL  =  ±0.03  ±  0.5 T aerL  21 . 

The  uncertainty  of  the  SODA  derived  cloud  or  aerosol  optical  thickness  is  similar  to  currently  well- 
established  remote  sensing  products  like  MODIS  and  could  in  theory  reach  the  accuracy  of  what  was  aimed 
by  the  Glory  mission  (±  0.02  on  optical  depth22).  The  advantage  of  SODA  is  as  that  the  error  sources  are 
totally  different  from  passive  remote  sensing,  it  can  provide  interesting  complementary  information. 

2.  Still  for  aerosols,  the  theoretical  error  bar  of  the  optical  depth  over  liquid  water  clouds  is  around  0.025.  This 
comes  from  the  high  theoretical  stability  of  the  lidar  ratio  in  liquid  water  clouds23.  A  recent  investigation24 
proposed  that  using  a  constant  liquid  water  clouds  lidar  ratio  in  the  algorithm  could  create  an  error  up  to  0.1 
in  the  optical  depth.  The  polarization  calibration  in  CALIPSO  version  3  seemed  to  create  an  artificial 
day/night  polarization  bias  of  around  30%25  which  would  impact  the  multiple  scattering  retrieval.  All  these 
effects  are  taken  into  account  and  mitigated  by  the  SODA  monthly  calibration  processes  done  prior  to  the 
optical  depth  retrieval.  Nevertheless,  because  of  the  novelty  of  the  method,  an  exact  error  budget  still  has  to 
be  derived  to  evaluate  how  far  is  the  data  product  from  the  theoretical  estimate.  Comparisons  with 
methodology  based  on  different  assumptions26, 27  will  provide  more  insight  in  the  matter. 

3.  For  cirrus  clouds  optical  depth  xcirL  over  the  ocean,  the  lack  of  well-established  reference  does  not  let  a  lot 
of  choice  but  to  be  cautious  and  provide  an  error  budget  that  is  likely  to  be  overestimated.  The  total  error 
we  have  determined  for  cirrus  optical  depth  based  on  remote  sensing  observations  is  equal  to 
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A  T drL  ~  (0.08)2  +  (0.25  TcirLf  (see  reference  [2]).  As  a  reference,  optical  depth  idoud  of  liquid  water 
clouds,  which  microphysical  properties  (spherical  water  droplets)  is  well  understood,  MODIS  uncertainty  is 
At doud  =  0.2t doud  •  However,  a  significantly  larger,  but  not  clearly  quantified,  uncertainty  exists  for  COD 

smaller  than  528, 29 .  Providing  an  uncertainty  estimate  for  the  optical  thickness  of  cirrus  clouds  derived  from 
passive  sensors  is  more  problematic  as  the  particles  forms,  shape  and  roughness  vary. 

4.  Over  the  ocean,  lidar  ratios  and  column  AOD  derived  from  CALIPSO  lidar  data  using  the  SODA  technique 
have  been  compared  with  those  measured  by  the  NASA  Langley  Research  Center  airborne  High  Spectral 
Resolution  Lidar  (HSRL)  for  several  cases  of  coincident  underflights19.  In  all  the  cases  studied  in  reference 
[19],  the  average  standard  deviation  between  the  SODA  and  HSRL  lidar  ratios  was  about  14%.  As  for 
liquid  water  clouds,  lidar  ratio  retrieval  from  space  is  a  novelty  and  more  work  would  be  useful  on  the 
topic.  Nevertheless,  SODA  is  as  far  as  we  know  the  only  global  scale  lidar  ratio  product  which  is  available 
to  the  scientific  community. 


4.  PRESENTATION  OF  THE  RESULTS 

Contrary  to  a  relatively  popular  belief,  performing  a  meaningful  comparison  with  AERONET  is  not  straightforward 
for  a  track  only  data  product  over  the  ocean.  Some  insight  of  the  difficulty  to  compare  a  track  only  product  with 
AERONET  is  provided  by  a  recent  study29  (for  both  land  and  ocean,  issues  are  magnified  over  the  ocean  due  to  the 
lower  number  of  data). 

We  presented  in  several  communications  the  comparison  between  SODA,  MODIS  and  POLDER.  Overall,  the 
agreement  is  good.  However,  the  intrinsic  differences  between  MODIS  and  POLDER  make  it  difficult  to  reach 
strong  conclusions.  In  general,  SODA  version  1  seemed  to  have  a  better  correlation  with  POLDER  than  with 
MODIS  (and  more  specifically  for  dust  particles)  as  determined  from  a  few  cases  but  the  reverse  is  found  in  the  data 
we  present  here  based  on  SODA  version  2  (Fig.  5).  Systematic  analyses  will  have  to  be  done  to  go  further.  For  the 
month  of  August  2010  (Fig.  5),  SODA  optical  depth  is  generally  slightly  higher  than  MODIS  and  lower  than 
POLDER.  It  may  be  due  to  the  presence  of  smoke  over  the  ocean  (Fig.  6)  which  absorption  properties  are  not 
always  easy  to  retrieve  by  passive  sensors.  Several  collocated  flight  with  the  NASA  LaRC  HSRL  system  were 
performed  during  this  time  period  and  supplemental  analysis  are  underway. 
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Fig.  5.  Upper  part:  MODIS  (collection  5)  optical  depth  as  a  function  of  SODA  optical  depth  (532  nm  left,  1064  nm 
right).  Bottom  part:  same  figure  using  the  POLDER  data.  All  aerosols  (cloud  filtered)  data  over  the  ocean  are  used. 

Figures  6,  7,  8  and  9  illustrate  some  results  that  were  obtained  from  SODA  version  2  by  filtering  the  data  with  the 
scene  flags.  The  optical  depth  of  dust  and  smoke  (as  classified  by  the  VFM)  are  shown  on  Fig.  6  over  the  ocean  and 
on  Fig.  7  over  liquid  water  clouds.  Fig.  8  and  9  show  the  same  data  but  for  the  lidar  ratio  retrieval.  Although  maps 
provide  an  easy  way  to  illustrate  the  data,  they  can  sometime  provide  a  false  sense  of  knowledge.  For  example,  the 
optical  depth  of  dust  over  the  liquid  water  clouds  is  very  low,  and  specific  attention  to  details  (or  alternatively  ensure 
that  there  is  a  significant  number  of  data  points)  would  be  required  before  reaching  conclusions  from  the  analysis  of 
the  optical  depth  and  lidar  ratio  in  that  case. 
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Fig.  6.  For  August  2010  nighttime  data.  532  nm  optical  depth  retrieval  for  dust  (left)  and  smoke  (right),  over  the 
ocean.  Data  have  been  filtered  with  the  SODA  scene  flags.  The  red  areas  represent  where,  based  on  their  location,  the 
features  appear  to  be  well  filtered  by  the  VFM. 
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Fig.  7.  Same  as  previous  figure  but  for  aerosols  above  liquid  water  clouds. 
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Fig.  8.  Same  as  Fig.  6  but  for  the  532  nm  lidar  ratio  of  dust  (left)  and  smoke  (right)  over  the  ocean. 


Fig.  9.  Same  as  Fig.  7  but  for  the  532  nm  lidar  ratio  of  dust  (left)  and  smoke  (right)  above  liquid  water  clouds. 


DISCUSSION  OF  THE  1064  NM  LIDAR  RATIO 

Before  concluding  this  paper,  we  would  like  to  provide  a  word  of  caution  about  the  1064  nm.  The  1064  nm  channel 
calibration  uncertainty  is  higher  than  the  532  nm.  Even  if  the  calibrations  processes  performed  at  the  different  stages 
of  the  algorithm  mitigate  the  different  error  sources,  there  are  less  established  references  at  this  wavelength 
(specifically,  it  is  more  difficult  to  validate  the  1064  nm  data  with  the  HSRL  underflights).  If  the  calibration 
processes  work  as  well  at  1064  nm  than  at  532  nm,  the  remaining  question  is  to  which  extent  the  baseline  of  optical 
depth  we  use15  is  accurate  at  1064  nm.  An  inaccurate  baseline  will  manifest  itself  as  a  bias  for  the  lowest  1064  nm 
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optical  depth.  If  this  bias  is  stable  (which  has  to  be  for  the  concept  of  baseline  to  be  valid),  it  will  be  identical  for 
different  atmospheric  features. 
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Fig.  10.  For  August  2010  nighttime.  1064  nm  optical  depth  as  a  function  of  the  532  nm  optical  depth  for  different  kind 
of  atmospheric  features.  Upper  part:  marine  (left)  and  dust  (right).  Bottom:  ice  (left)  and  liquid  (right)  water  clouds. 

Here  we  illustrated  the  concept  with  coarse  mode  aerosols  and  clouds  (Fig.  10)  as  the  analysis  is  more 
straightforward  when  we  expect  some  level  of  wavelength  independency.  SODA  allows  to  derive  the  effective  OD 
of  semi-transparent  layers  but  the  multiple  scattering  contribution  has  to  be  further  corrected  (note  that  no  multiple 
scattering  corrections  are  included  in  SODA).  For  large  particles  like  liquid  water  droplets  and  ice  cloud  crystals,  we 
can  see  the  effect  of  the  forward  peak  and  of  the  differential  geometry  of  the  lidar  emission/reception.  The  laser  light 
is  within  the  telescope  field  of  view  at  532  nm  but  the  laser  divergence  is  larger  than  the  telescope  divergence  at 
1064  nm  (see  CALIPSO  level  1  ATBD  table  2.1  note  a).  When  this  effect  is  combined  with  the  larger  size  of  the 
particle  forward  peak  at  1064  nm,  more  light  escape  the  telescope  at  1064  nm  than  at  532  nm  which  creates  a 
slightly  larger  attenuation  (and  optical  depth).  This  effect  is  not  noticeable  for  smaller  aerosol  particles  like  dust  and 
marine  aerosols.  In  version  3,  smoke  aerosols  were  sometime  classified  as  marine  aerosols.  It  will  typically  create 
the  “two  branches”  visible  on  Fig.  10  for  marine  aerosols. 

Although  we  do  not  want  to  provide  a  definitive  number  from  this  preliminary  analysis  of  one  month  of  data,  but  it 
is  quite  conceivable  that  the  1064  nm  optical  depth  are  slightly  underestimated  (about  10%)  with  respect  to  532  nm 
ones.  532  nm  seems  slightly  larger  at  the  lowest  optical  depth  values  on  Fig.  10.  For  users  interested  in  the  1064  nm 
lidar  ratio  of  optically  thin  features,  we  suggest  a  similar  analysis  to  be  performed  on  the  optical  depth  to  improve 
the  uncertainty  estimate. 


CONCLUSION  AND  PERSPECTIVES 

We  hope  that  this  short  document  will  help  on  enlight  the  reader  on  the  features  of  the  SODA  product,  its 
advantages  and  current  limitations. 
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In  term  of  short  term  perspectives,  further  analysis  of  the  atmospheric  product  (polarization  and  intensity)  will 
provide  insight  on  atmospheric  fine  mode  fraction  and  on  the  use  of  CloudSat  as  a  reference  to  correct  the 
underwater  lidar  polarization  observations  which  has  different  applications  to  ocean  science6. 

Long  term  perspectives  are  more  difficult  to  foresee.  If  we  look  at  the  past  few  years,  SODA  went  from  what  was  a 
concept  designed  to  provide  the  aerosol  optical  depth8  to  applications  to  aerosols  and  clouds  optical  properties 
(optical  depth,  color  ratio  and  lidar  ratio),  water  vapor  absorption  for  W-band  radar  and  marine  biology. 
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